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ABSTRACT 



We present more than three years of observations at different frequencies, from radio to high-energy y-rays, of the Narrow-Line 
Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z = 0.585). This source is the first NLS1 detected at energies above 100 MeV and 
therefore can be considered the prototype of this emerging new class of y-ray emitting active galactic nuclei (AGN). The observations 
performed from 2008 August 1 to 201 1 December 31 confirmed that PMN J0948+0022 generates a powerful relativistic jet, able to 
develop an isotropic luminosity at y-rays of the order of 10 48 erg s" 1 , at the level of powerful quasars. The evolution of the radiation 
emission of this source in 2009 and 2010 followed the canonical expectations of relativistic jets, with correlated multiwavelength 
variability (y-rays followed by radio emission after a few months), but it was difficult to retrieve a similar pattern in the light curves 
of 2011. The comparison of y-ray spectra before and including 2011 data suggested that there was a softening of the high-energy 
spectral slope. We selected five specific epochs to be studied by modelling the broad-band spectrum, characterised by an outburst at 
y-rays or very low/high flux at other wavelengths. The observed variability can largely be explained either by changes in the injected 
power, the bulk Lorentz factor of the jet or the electron spectrum. The characteristic time scale of doubling/halving flux ranges from 
a few days to a few months, depending on the frequency and the sampling rate. The shortest doubling time scale at y-rays is 2.3 ± 0.5 
days. These small values underline the need of highly-sampled multiwavelength campaigns to better understand the physics of these 
sources. 

Key words, galaxies: individual: PMN J0948+0022 - galaxies: Seyfert - galaxies: jets - gamma-rays: galaxies 



1. Introduction 

PMN J0948+0022 (also known as SDSS J094857.3 1 +002225.4, 
Z = 0.585 ± 0.001Q) is a very peculiar AGN, with an optical 
spectrum of a Narrow-Line Seyfert 1 (NLS1) galaxy (Williams 

1 Sloan Digital Sky Survey Team, SDSS Data Release 2, Nov. 18, 
2004 from: http://www.sdss.org/dr2/products/spectra/getspectra.html 



et al. 2002). Zhou et al. (2003) first reported the unexpected con- 
nection of its NLS 1 character with strong and fiat-spectrum ra- 
dio emission. The optical spectrum of PMN J0948+0022 ex- 
hibits the usual characteristics of NLS Is (cf. Osterbrock & 
Pogge 1985, Goodrich 1989): FWHM(HyS) = 1500 ± 55 km/s, 
[OIII]/H/3 < 3, strong optical Fell emission (Zhou et al. 2003). 
Such types of AGN are generally radio-quiet (Ulvestadt et al. 
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1995, Komossa et al. 2006), so it was a surprise to detect strong 
radio emission. The radio power at 5 GHz was measured to be 
~ 10 43 erg s , a relatively small value for a blazar, but very high 
for a NLS1, resulting in a radio loudness RL = fs ghz//440 nm ^ 
10 3 (Zhou et al. 2003). This value is among the highest for a 
NLS1 (although not the highest, cf. Fig. 3 in Foschini 2011a). 

The inverted spectrum (as-io ghz 0.24, with f v cc v~ a ) and 

high brightness temperature (7b > 10 13 K) suggest the presence 
of a relativistic jet viewed at small angles (Zhou et al. 2003, Doi 
et al. 2006), similar to the conditions found in blazars (Blandford 
& Rees 1978). 

Confirmation of a powerful relativistic jet viewed at small 
angle arrived with the launch of the Fermi Gamma-ray Space 
Telescope (hereafter Fermi) carrying onboard the Large Area 
Telescope (LAT, LAT Coll. 2009a), which operates in the en- 
ergy band from 100 MeV to more than 300 GeV and performs 
an all-sky survey every two orbits (~ 3 hours). During its early 
months of operation, LAT detected variable high-energy y-ray 
radiation from PMN J0948+0022 (LAT Coll. 2009b,c, 2010), 
and after one year the number of y-NLSls increased to four 
(LAT Coll. 2009d). To date there are five high-confidence (Test 
Statistic TS > 25, see Mattox et al. 1996 for a definition of 
TS) and three low-confidence (TS = 9- 25) y-ray detections of 
NLS 1 plus a similar number of candidates in the possible parent 
population^ (see Foschini 201 la for a review). 

The importance of the research resides in understanding a 
few fundamental points. On one side, there are blazars (beamed 
population) plus radio galaxies (unbeamed parent population): 
they have the mass of the compact object in the range ~ 10 8 - 
10 10 M G , accretion discs spanning a wide range of luminosities, 
from ~ 10~ 6 to a significant fraction of the Eddington limit (e.g. 
Ghisellini et al. 2010, Foschini 2011b), and the host galaxy is 
elliptical (e.g. Kirhakos et al. 1999). At the other side, there are 
the y-NLS 1 s (beamed population) plus an unbeamed parent pop- 
ulation yet to be exploited, which have the central supermassive 
black hole with masses between ~ 10 6 - 1O 8 M , accretion disc 
luminosities close to the Eddington limit (e.g. Foschini 201 la,b), 
and a few of them are hosted in spirals possibly having under- 
gone a recent merger (see Yuan et al. 2008, Foschini 2011a,b, 
2012, Hamilton & Foschini 2012). Therefore, it is now possible 
to study AGN with powerful relativistic jets in objects that ex- 
hibit an unexplored range of central black hole masses and rates 
(cf Fig. 8 in Foschini 201 la). 

PMN J0948+0022 has been highly active during the past 
three and a half years. Radio data from Owens Valley Radio 
Observatory (OVRO) extending back to one year before the 
launch of Fermi (see LAT Coll. 2009b, Fig. 2, Panel C) displayed 
a radio flux density at 15 GHz that reached up to ~ 1 Jy on 
2008 January, a factor ~ 2 - 3 above the value at the time of the 
discovery (2008 August-December). The first multiwavelength 
(MW) campaign organised in 2009 March- July revealed contin- 
uing activity (LAT Coll. 2009c, see also Giroletti et al. 201 1 for 
e-VLBI observations) and in 2010 July the source erupted into 
an exceptional outburst with a peak y-ray isotropic luminosity 
of ~ 10 48 erg s 1 (LAT Coll. 2010b, Foschini 2010, Foschini et 
al. 201 la,b). The new data presented here show that the activity 
continued into 201 1 as well. 

PMN J0948+0022 has been the target of several observations 
at almost every waveband during the past three and a half years 
(2008 August 1 - 201 1 December 31), resulting in an impressive 



2 The current status of the search for high-energy y-ray emission 
from NLS Is and the related suggested parent population is available 
at the web page http: //tinyurl . com/gnlsls prepared by LF. 



wealth of data. We present here a comprehensive study of all 
these data in order to set up some firm pillars in our understand- 
ing of the prototype of this new class of y-ray AGN. Part of the 
data presented here has been already published by the LAT Coll. 
(2009b,c,d), Foschini et al. (2011a,b), Angelakis et al. (2012), 
but they have been reanalysed by using the most recent software 
version and calibration database. This is particularly important 
in the case of LAT data, where the change from the P6 to P7 in- 
strument response function determined important changes in the 
analyses (cf. LAT Coll. 2010c with LAT Coll. 2012). The effects 
of these changes on the 2010 July outburst have been presented 
in Foschini et al. (201 lb), and the present work reports about the 
reanalysis of the whole period 2008-201 1. In addition, the 201 1 
data are presented here for the first time. 

The work is structured as follows: after the introduction, 
Sect. 2 deals with the descriptions of the individual instruments 
and facilities used for the observations, and the techniques for 
data analysis. A general overview of the data (avoiding interpre- 
tations) is given in Sect. 3. Five specific periods have been se- 
lected and the corresponding broad-band spectra have been stud- 
ied with the one-zone synchrotron and inverse-Compton model 
by Ghisellini & Tavecchio (2009). This part is described in 
Sect. 4. Some final remarks (Sect. 5) conclude the work. 

Throughout this work, the luminosities have been calculated 
by assuming a ACDM cosmology with a Hubble-Lemaitre con- 
stant H = 70 km s -1 Mpc -1 and Q A = 0.73 (Komatsu et al. 
2011). The luminosity distance of PMN J0948+0022 is then 
3462 Mpc and 1 arcsec corresponds to 6.68 kpc. In the fol- 
lowing, the dates will be indicated either with calendar scheme, 
or Modified Julian Date (MJD), or days since 2008 August 1 
(MJD-54679). 

2. Data Analysis 

2.1. Fermi Large Area Telescope 

LAT photons of class 2 ("source" type) with energy between 
100 MeV and 100 GeV recorded between 2008 August 4 15:43 
UTC and 2011 December 31 24:00 UTC in a circular region 
around the position of PMN J0948+0022 with 10° radius were 
retrieved from the Fermi Science Support Centefl The data were 
analysed by using LAT Science Tools v. 9 . 27 . 1, together 
with the Instrument Response Function (IRE) Pass 7 and the cor- 
responding isotropic and Galactic diffuse background modelfl 
The adopted procedures are standard and described in detail in 
the analysis threads available online at the Fermi SSC. 

The data were fitted with a single power-law model. The 
analysis of all the data integrated into one single block resulted 
in 4951 photons predicted by the model, which in turn corre- 
sponds to a flux Fo.i-ioo GeV = (1.36 ± 0.03) x 10~ 7 ph crrr 2 s _1 
and a photon index F = 2.67 ± 0.03. The likelihood of detec- 
tion is expressed by means of the Test Statistic (TS , Mattox et 
al. 1996), which, as a rule of thumb, is linked to the significance 
no- as \[f~S cr. For the Fermi observations of PMN J0948+0022, 
TS =2015. 

In the 2FGL catalogue (LAT Coll. 2012), the energy distri- 
bution N(E) of PMN J0948+0022 is fitted with a log-parabola 
model of the form: 

N(E) = No^- j (1) 

3 http : //f ermi .gsfc.nasa.gov/ssc/data/ 

4 See the documentation available online at the Fermi SSC. 
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Fig. 1. y-ray power-law model photon index as a function of 
the 0.1 - 100 GeV flux. The data points represent 7-day bins. 
The dashed line marks the FermifLPZY 7-day detection threshold 
(TS = 9) as extrapolated from LAT Coll. (2010c). 



where A^o is the normalisation [ph cm 2 s 1 MeV ], or is the 
photon index at the pivot energy E p [MeV], and /3 is the curvature 
index. Specifically for PMN J0948+0022, No = (1.39 ± 0.07) x 
10~ 10 ph cm- 2 s- 1 MeV- 1 , a = 2.26 ± 0.08,/3 = 0.26 + 0.06, and 
E p ~ 272 MeV and refer to the period 2008 August 4-2010 
July 31 (LAT Coll. 2012). 

For comparison, using the same model for the integrated data 
considered here (2008 August 4-2011 December 31) we ob- 
tain values of N Q = (1.23 ± 0.03) x 10" 10 ph cirT 2 s _1 MeV" 1 , 
a = 2.45+0.03,/? = 0.23+0.02, and E p = (3 13 ±3) MeV. In total, 
4802 photons are predicted by the model, corresponding to an in- 
tegrated flux of F . i-ioo Gev = (1.21 + 0.03) x 10~ 7 ph cm" 2 s"\ 
and TS = 2070. The changes are basically due to the addition 
of the 201 1 data, when the source was strongly active for almost 
the entire year (see Section 3, Fig. [3] top panel), and slight soft- 
ening of the spectrum occurred (a increased by ~ 0.2, while the 
curvature /3 remained constant). 

Fig-Edisplays the behaviour of the y-ray spectral slope as a 
function of the 7 -day integrated flux in the 0. 1 - 100 GeV energy 
band. The measured photon index is that of the power-law model 
and no strong spectral changes can be observed on short time 
scales, while there could be on longer periods, as shown from 
the above comparison of the 2FGL data with the present work. 

The y-ray light curve in Fig.[3](fo/? panel) and I A. II (top left 
panel) was built with one-day time bins, assuming a power-law 
model with T = 2.67 (i.e. the value from the integration of the 
whole data set), and with a detection threshold of TS = 9. Since 
the source is a confirmed y-ray emitter, we relaxed the statis- 
tical requirements for valid detections that has been adopted in 
recent works (e.g. Foschini et al. 2011a,b). Due to its scanning 
mode, Fermi monitors any source in the sky almost continuously. 
Sometimes, PMN J0948+0022 was not detected on daily time- 
scale, but for clarity we do not plot upper limits. The TS — 4 av- 
erage upper limit for one-day exposure (~ 2.3x 10~ 7 ph cirT 2 s _1 , 
by fixing F = 2.67) is indicated in Fig.[3](fo/? panel) and lA.l \ (top 
left panel) with an horizontal dotted line. 



2.2. Swift 

The data of PMN J0948+0022 from the three Swift instruments, 
the Burst Alert Telescope (BAT, Barthelmy et al. 2005), the X- 



Ray Telescope (XRT, Burrows et al. 2005), and the Ultra Violet 
Optical Telescope (UVOT, Roming et al. 2005) were down- 
loaded from the HEASARC archive and analysed with the Swift 
tools included in the HEASoft v. 6 . 12 packag^l together with 
the CALDB updated on 20 1 2 March 22 . 

Swift observed the source for the first time on 2008 
December 5 soon after the discovery of the y-ray emission (LAT 
Coll. 2009b). During the 2009 MW campaign, there were 1 1 
snapshots (LAT Coll. 2009c) and one a few days before the out- 
burst of 2010 July (Foschini et al. 2011a). Ten more exposures 
were performed in 201 1, as part of a monitoring program linked 
to the Effelsberg radio observations. In total, from the period 
2008-201 1, there are 23 observations with XRT exposures from 
1 to 5 ks each. 

All the available pointed BAT observations were combined 
to search for a possible detection. In total, ~ 86 ks of data 
were accumulated, but no detection was made. The upper lim- 
its at 99.97% confidence level were 1.7 x 10~ 10 erg cirT 2 s _1 
in the 20 - 40 keV energy band and 2.3 x 10~ 1() erg cirT 2 s _1 
in the 40 - 100 keV band. Cusumano et al. (2010) reported a 
detection by integrating 54 months of BAT data (thus including 
also the Swift observations when PMN J0948+0022 was not on 
axis, but is still within the BAT field of view), with a flux of 
(1.2 + 0.8) x lO" 11 erg cirT 2 s" 1 in the 10 - 150 keV energy 
bancfl Despite the large error, the flux seems to be consistent 
with what is expected from the modelling of the spectral energy 
distribution (see Section 4, Fig. [4). 

XRT was set to work in photon counting mode (Hill et al. 

2004) and we analysed all the single-to-quadruple events (grades 
0-12). The extracted spectra were rebinned to have at least 20 
counts per bin to apply the x 2 statistics. When this was not pos- 
sible, we evaluated the likelihood using Cash statistics (Cash 
1979). The data were fitted with a redshifted power-law model 
with Galactic absorption (Nn = 5.23 x 10 20 cirT 2 , Kalberla et al. 

2005) . The individual snapshot does not show significant spec- 
tral changes as a function of the emitted power, but the statistics 
are not sufficient to test if more complex models can be em- 
ployed (Fig. O. When all the X-ray data are integrated (the re- 
sulting global exposure is 81.6 ks), the spectrum is best fitted 
(X 2 = 241.66, dof = 200) with r = 1.67 + 0.03 and normalisa- 
tion at 1 keV equal to (1 .50+0.05)x 10~ 3 ph cirT 2 s 1 keV 1 . The 
integrated observed flux in the 0.3 - 10 keV band is (4.5 + 0.2) x 
10~ 12 erg cm" 2 s _1 . The modelling of the integrated spectrum 
can be improved (x 2 = 210.54, dof = 198, F-test > 99.99%) by 
adopting a broken power-law model with the following param 
eters: V x = 1.91 + 0.10, T 2 = 1.56!°^, £ break 



1.22^,* keV, 
1 at 1 keV. 



and normalisation (6.7 + 0.3) x 10 ph cm s keV 

It is possible to model the spectrum also by using log- 
parabolic models (logpar in the xspec syntax), which have 
been applied to blazars, specifically to High-Frequency-Peaked 
BL Lac Objects (e.g. Massaro et al. 2004, Foschini et al. 2007). 
The best-fit parameters are a = 1.79 + 0.04 and/? = -0.27 + 0.04, 
with the pivot energy fixed at 1 keV. The x 2 test resulted in a 
value of 210.34 for 199 degrees of freedom, which in turn trans- 
lates into a F-test better than 99.99% with respect to the single 
power-law model. A variant of the above model, eplogpar in 
xspec (see Tramacere et al. 2007), resulted in similar values for 
the^ 2 and F- tests (x 2 = 210.46, dof = 199, F-test > 99.99%) 



5 Including the XRT Data Analysis Software (XRTDAS) developed 
under the responsibility of the ASI Science Data Center (ASDC), Italy. 

6 It is worth noting that the 58-month survey by 
Baumgartner et al. does not report this detection. See: 
http : //swift . gsfc .nasa . gov/docs/swift/results/bs58mon/ 
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Fig. 2. X-ray power-law model photon index plotted as a func- 
tion of the 0.3-10 keV luminosity. No obvious luminosity depen- 
dent spectral changes are seen. 



and the following model parameters: the peak energy in vF v is 
2.4 ± 0.2 keV and the curvature B = -0.28 + 0.09. In both cases, 
the improvement with respect to the single power-law model is 
significant, but there is no preference with respect to the broken 
power-law model. Basically, the spectrum indicates a change in 
the slope around 1-2 keV, with the high-energy photon index 
harder than that at lower energies. 

UVOT data were taken with all six filters (V, B, U, UVWl, 
UVM2, UVW2) for most the observations. They were analysed 
by extracting the source counts from a 5"-sized region and the 
background from a source-free annulus centred on the source, 
and with internal and external dimensions of 7" and 40", re- 
spectively. The observed magnitudes were dereddened by using 
A v = 0.28 and the extinction laws by Cardelli et al. (1989). The 
intrinsic magnitudes were then converted in physical units by 
using the zero-points and conversion factors of the Swift CALDB 
(Poole et al. 2008, Breeveld et al. 2010). 

2.3. Instituto Nacional de Astrofisica, Optica y Electrdnica 
(INAOE) 

PMN J0948+0022 was observed at near-infrared (NIR) wave- 
lengths at the 2.1 m telescope "Guillermo Haro" (Cananea, 
Sonora, Mexico). The NIR camera "CANICA" has a Rockwell 
1024 x 1024 pixel Hawaii infrared array (0.32 arcsec/pix), op- 
erating at 75.4 K, with standard J (1.164 - 1.328 //m), H 
(1.485 - 1.781 fim) and K s (1.944 - 2.294 jum) filters in place. 
Observations were carried out in series of 10 dithered frames in 
each filter, with a proper number of additional observations for 
the K s filter. Data sets were corrected for bias and flat-fielding, 
the latter obtained from sky frames derived from the dithered 
ones, and coadded. 

2.4. Effelsberg and IRAM 

The cm/mm radio light curves of PMN J0948+0022 have been 
obtained within the framework of a Fermi related monitoring 
program of y-ray blazars (F-GAMMA program, Fuhrmann et 
al. 2007, Angelakis et al. 2008) as well as a dedicated NLS1 
monitoring program. The millimetre observations are closely co- 
ordinated with the more general flux density monitoring con- 
ducted by IRAM (Instituto de Radioastronomia Milimetrica), 



and data from both programs are included in this paper. The 
overall frequency range spans from 2.64 GHz to 142 GHz us- 
ing the Effelsberg 100-m and IRAM 30-m (at Pico Veleta) tele- 
scopes. 

The Effelsberg measurements were conducted with the sec- 
ondary focus heterodyne receivers at 2.64, 4.85, 8.35, 10.45, 
14.60, 23.05, 32.00, and 43.00GHz. The observations were per- 
formed quasi-simultaneously with cross-scans, that is slewing 
over the source position, in azimuth and elevation direction with 
adaptive number of sub-scans for reaching the desired sensitiv- 
ity (for details, see Fuhrmann et al. 2008; Angelakis et al. 2008). 
Consequently, pointing off-set correction, gain correction, atmo- 
spheric opacity correction and sensitivity correction have been 
applied to the data. 

The IRAM 30-m observations were carried out with cali- 
brated cross-scans using the ABCD' SIS (until March 2009) 
and new EMIR horizontal and vertical polarisation receivers op- 
erating at 86.2 and 142.3 GHz. The opacity corrected intensities 
were converted into the standard temperature scale and finally 
corrected for small remaining pointing offsets and systematic 
gain-elevation effects. The conversion to the standard flux den- 
sity scale was done using the instantaneous conversion factors 
derived from frequently observed primary (Mars, Uranus) and 
secondary (W3(OH), K3-50A, NGC7027) calibrators. 

The complete data set has been published in Angelakis et 
al. (2012), together with data of other NLSls. The brightness 
temperature of PMN J0948+0022 was measured as 8 x 10 12 K 
at 4.85 GHz, 2 x 10 12 K at 14.6 GHz, and 1.5 x 10 11 K at 32 
GHz, which can be translated into requirement of Doppler fac- 
tors greater than 6, 4, and 2, respectively (Angelakis et al. 2012). 

2.5. Metsahovi 

The Metsahovi radio telescope is a radome enclosed paraboloid 
antenna (13.7 m) situated in Finland and operating at 37 GHz. 
The detection limit (signal-to-noise ratio > 4) of the telescope at 
37 GHz is of the order of 0.2 Jy under optimal conditions and 
with a typical integration time of 1200- 1400 s. The flux density 
scale was set by observations of DR 21, while 3C 84, 3C 274, 
and NGC 7027 were used as secondary calibrators. More details 
about data reduction and analysis can be found in Terasranta et 
al. (1998). The error in the flux density takes into account the 
background and the absolute calibration. 

2.6. Owens Valley Radio Observatory 

Regular 15 GHz observations of PMN J0948+0022 were carried 
out as part of a high-cadence y-ray blazar monitoring program 
using the OVRO 40 m telescope (Richards et al. 201 1). This pro- 
gram, which commenced in late 2007, now includes about 1600 
sources, each observed with a nominal twice per week cadence. 

The OVRO 40 m uses off-axis dual-beam optics and a cryo- 
genic high electron mobility transistor (HEMT) low-noise am- 
plifier with a 15.0 GHz center frequency and 3 GHz bandwidth. 
The total system noise temperature is about 52 K, including re- 
ceiver, atmosphere, ground, and CMB contributions. The two 
sky beams are Dicke switched using the off-source beam as a ref- 
erence, and the source is alternated between the two beams in an 
ON-ON fashion to remove atmospheric and ground contamina- 
tion. A noise level of approximately 3-4 mJy in quadrature with 
about 2% additional uncertainty, mostly due to pointing errors, 
is achieved in a 70 s integration period. Calibration is achieved 
using a temperature-stable diode noise source to remove receiver 
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2009 Jan 1 May 5-15 2010 Jan 1 Jul 8 201 1 Jan 1 Jun 21 Oct 1 1 




Days since 2008 Aug 1 [MJD - 54679] 

Fig. 3. Light curves at various frequencies. From top to bottom: y-rays 0.1 - 100 GeV from FermifLPC[ with 1 day time bin 
[10- 7 ph cm 2 s" 1 ]; X-ray 0.3 - 10 keV from Swift/XRT [c s" 1 ]; U filter (350 nm) from Swift/UVOT [mJy]; H filter (1.65 /mi) from 
INAOE [mJy]; 86 GHz from IRAM [Jy]; 37 GHz from Metsiihovi [Jy]; 15 GHz from OVRO, Effelsberg, and MOJAVE [Jy]; 5 GHz 
from Effelsberg and Medicina. Time starts on 2008 August 1 00:00 UTC (MJD 54679). Fermi performs an all-sky survey every two 
orbits (3 hours) and therefore it observed PMN J0948+0022 almost continuously. However, for clarity, we do not plot the upper 
limits. The TS — 4 average sensitivity of LAT for one day exposure is indicated with the horizontal dotted line in the top panel, 
which corresponds to ~ 2.3 x 10~ 7 ph cirT 2 s _1 (assuming T = 2.67). Vertical dashed lines are shown at the beginning of each year 
(Jan 1). Vertical dot-dashed lines indicated the days referring to the spectral energy distributions selected as example (see Sect. 4). 
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gain drifts and the flux density scale is derived from observations 
of 3C 286 assuming the Baars et al. (1977) value of 3.44 Jy at 
15.0 GHz. The systematic uncertainty of about 5% in the flux 
density scale is not included in the error bars. Complete details 
of the reduction and calibration procedure are found in Richards 
etal. (2011). 

2. 7. VLBA 2 cm (MOJAVE Program) 

The MOJAVE Program (Monitoring Of Jets in Active galac- 
tic nuclei with VLBA Experiments, Lister et al. 2009) is 
a high-resolution high-dynamic range (about 8000:1) survey 
at 15.4 GHz performed with the Very Large Baseline Array 
(VLBA). It counts about 300 active galactic nuclei in the 
Northern hemisphere, which are observed on a roughly regu- 
lar basis. Presently, the program includes four radio-loud NLS 1 s 
detected at y-rays (J0324+3410, J0849+5108, J0948+0022, 
J1505+0326). Specifically, PMN J0948+0022 has been ob- 
served 11 times between 2009 May 28, and 2011 December 12 
(Table[TJ. Details of data processing are available in Lister et al. 
(2009) and Lister & Homan (2005). The analysis of the 2009 
and part of 2010 data resulted in a brightness temperature of the 
core in excess of 6 x 10 12 K, a jet apparent opening angle of 21°, 
a jet position angle of 24°, an average fractional polarisation of 
0.8%, and an average Electric Vector Position Angle (EVPA) of 
142° (Lister et al. 201 1). There is negligible difference between 
the single dish (OVRO) and MOJAVE (VLBA) 2 cm flux densi- 
ties at all epochs, indicating no bright arcsec scale radio structure 
(Lister etal. 2011). 

The three initial epochal showed two moving features, but 
these faded rapidly and thus their measured speeds cannot be 
considered reliable. 

Table 1. VLBA observations performed within the MOJAVE 
project. 



Date 


Intensity [mJy] 


Polarisation [%] 


EVPA [i 


2009 May 28 


384 


0.7 


131 


2009 Jul 23 


462 


1.2 


146 


2009 Dec 10 


675 


0.5 




2010 Sep 17 


669 


3.7 


51 


2010 Nov 04 


627 


1.8 


49 


2010 Nov 29 


486 


1.6 


49 


2011 Feb 20 


498 


1.0 


19 


2011 May 26 


540 


1.8 


49 


2011 Jun24 


389 


2.1 


51 


2011 Sep 12 


343 


1.3 


40 


2011 Dec 12 


439 


0.8 


77 



2.8. Medici na 

We observed PMN 0948+0022 with the 32m Medicina radio 
telescope eight times between 2011 June and November, at a 
frequency of 5 and/or 8.4 GHz. At each epoch and frequency, 
we performed ~ 10 cross scans on the source using the new 
Enhanced Single-dish Control System (ESCS) acquisition sys- 
tem. The typical on source integrated time is ~ 1 minute per 
band. Since the signal-to-noise ratio in each scan across the 

7 Images with animations can be found at: 
|http : //tinyurl . com/mo j aveS948pSQ22 



source was low (typically ~ 3), we performed a stacking anal- 
ysis of the scans, which allowed us to significantly improve 
the signal-to-noise ratio and the accuracy of the measurement. 
Finally, we calibrated the flux density with respect to simultane- 
ous observations of 3C 286. 

2.9. Sloan Digital Sky Survey (SDSS) 

We also reanalysed the non-contemporaneous data of the SDSS. 
The images were taken in 1999 March, and the spectrum (cover- 
ing 3800 - 9200 A, or rest-frame 2400 - 5600 A) in 2000 March. 
The host galaxy is unresolved in the SDSS images, which is not 
surprising, given its redshift (z = 0.5846). The spectrum shows 
evidence for a young stellar population, in addition to the nu- 
clear source (Hamilton & Foschini 2012). High-resolution ob- 
servations (e.g. with Hubble Space Telescope) are required to 
study the morphology of the host galaxy. 

3. Overview of the data 

All the data available for the present study are displayed in the 
figures of Appendix A, while a sub-sample with the most rep- 
resentative light curves is shown in Fig. [3] The period covered 
is from 2008 August 1 to 2011 December 31. Both samples 
clearly show a high degree of activity in PMN J0948+0022 at 
all wavelengths. The most complete coverage is at y-rays with 
Fermi/LAJL (daily), and at 15 GHz with OVRO (twice per week), 
plus contributions from Effelsberg (once per month) and the 
MOJAVE Program (Table [TJ- At the other wavelengths the data 
are basically clustered around two main MW campaigns organ- 
ised in 2009 (LAT Coll. 2009c) and 201 1 (new data). 

3.1. 2008 

As already published in LAT Coll. (2009b), radio observations 
at OVRO (15 GHz) before the launch of Fermi indicated a strong 
radio emission with a peak of 0.81 Jy on 2008 Jan 9 (MJD 
54474). At the start of the scientific operations of the LAT on- 
board Fermi, PMN J0948+0022 was detected at a level of a few 
xl0~ 7 ph cm~ 2 s _1 making it the first NLS1 detected at high- 
energy y-rays (LAT Coll. 2009b, 2010a). 

3.2. 2009 

We confirm the findings of the 2009 MW Campaign, performed 
between 2009 Mar 20 and Jul 5 (LAT Coll. 2009c). The source 
displayed some activity at y-rays in April, with the most sig- 
nificant (i.e. with greatest TS) peak detected on April 1 (MJD 
54922; day 243 in Fig. [3) with a flux of (4.2 ± 1.5) x 10~ 7 
ph crrT 2 s _1 in the 0.1 - 100 GeV energy band with a photon in- 
dexT = 2.1+Q.2(TS = 43). Other detections on daily basis were 
recorded by ferm//LAT during almost the whole month, with 
fluxes of a few times 10 -7 ph crrr 2 s _1 (0.1 - 100 GeV). Near- 
Infrared observations at INAOE show a relatively bright state, 
with observed magnitudes / = 15.87 + 0.05, H = 14.90 + 0.07, 
K s = 13.99 + 0.09, while optical-polarimetric observations 
(V filter) performed in the period 2009 Mar 30 - Apr 10 at 
the KANATA telescope resulted in a 19% polarisation fraction 
(Ikejiri etal. 2011). 

Radio observations in 2009 April reported low flux densities 
of the order of ~ 0.2 Jy at 15 GHz (see Fig. [3}, but increasing 
and peaking in mid-end May. At 15 GHz, with the densest sam- 
pling, the peak of 0.43 Jy was measured on 2009 May 21 (MJD 
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54972, day 293 in figures), that is 50 days after the y-ray peak. 
Metsahovi (37 GHz) reported 0.83 Jy on the evening of 2009 
May 17 (MJD 54968.72, day ~ 290). Instead, on 2009 May 27 
(MJD 54978.77 or day ~ 300), Effelsberg measured flux den- 
sities reach 0.269 Jy at 10.45 GHz), and 0.189 Jy at 8.35 GHz 

(Fig.E2). 

The Swift optical-to-X-ray monitoring was much coarsely 
sampled. However, two observations performed on 2009 May 
5 and 15 show a drop in the 0.3-10 keV flux from (7.5 ± 0.3) x 
10~ 12 erg cirT 2 s to (2.1+0.4)xl0~ 12 erg crrT 2 s , correspond- 
ing to an intrinsic (i.e. redshift corrected) halving time scale of 
3.4 + 0.7 day. At the same time, the UVOT spectrum exhibited 
a slope change. Over a 5 -day period centred on 2009 May 5, 
when the y-ray flux was at level of (2.2 ± 0.4) x 10 -7 ph crrT 2 s _1 
(TS = 10), the optical/UV spectrum was rather flat. Ten days 
later the y-ray flux fell below the LAT detectability (< 8 x 10" 8 
ph cirT 2 s _1 , upper limit at TS =4 with 5 days integration). At 
the same time, a flux density drop of ~ 2.8 was seen in opti- 
cal filters, whereas the flux density decreased by only ~ 1.8 in 
the UV. Interestingly, the H filter detected a flux density drop on 

2009 May 1 (MJD 54952, day 273), prior to the optical-to-X-ray 
fluxes. The NIR flux density had decreased by roughly a factor 
2 since the measurements performed in 2009 April. 

The source ceased detectable y-ray activity from 2009 
August to about 2010 June, with the exception of a couple of 
detections in 2009 mid December. There was no information at 
optical-to-X-rays, but the flux density at radio frequencies was 
high as measured at 86 GHz (IRAM, more than 1 Jy on 2009 Dec 
7) and at 15 GHz (0.46 Jy on 2009 Dec 31). Interestingly, the 
VLBA 2-cm observation performed on 2009 Dec 10 (MOJAVE), 
resulted in the lowest polarisation fraction ever measured (0.5%, 
Table[T|). It was not possible to measure the EVPA in that epoch, 
but the following observation performed on 2010 Sept 17, re- 
vealed that it had changed significantly, rotating by about 90°, 
from 146° on 2009 Jul 23, to 51° (TableQ). 

3.3. 2010 

Then, on 2010 July 8 (day 706), PMN J0948+0022 exploded 
in its first known y-ray outburst, with isotropic luminosity at y- 
rays in excess of 10 48 erg s , comparable with powerful quasars 
(LAT Coll. 2010b, Foschini 2010, Foschini et al. 2011a,b). At 
that time, it was not possible to set up an efficient campaign, 
because of the close apparent position of the Sun, but the data 
available show a relatively high level of optical-to-X-ray fluxes. 
Measurements at radio frequencies show an increasing flux den- 
sity that reached the peak at 15 GHz about 46 days after the first 
y-ray outburst (0.55 Jy on 2010 Aug 22, day 752) and remain- 
ing almost constant at this level for about one hundred days. On 

2010 Sept 17 (day 777), the VLBA 2-cm observation reported 
the highest fractional polarisation measured to date at this wave- 
length (3.7%, Table [B- The flux density at 37 GHz had a local 
maximum on June 27 (0.88 Jy, day 695) and July 10 (0.81 Jy, 
day 708). The next measurement was only on Sept 2 with 0.5 Jy 
(day 762). The 86 GHz flux density was again at ~ 1 Jy level on 
Aug 2 (day 731). 

3.4. 2011 

The year 2011 data are characterised by a prolonged activity 
both at y-rays and radio frequencies. Specifically, Effelsberg and 
Medicina reported a 5 GHz flux density almost double the val- 
ues seen in 2008-2009. Starting from 2011 Feb 1 (day 914), 



Metsahovi (37 GHz) registered a flux density increase, reach- 
ing the maximum value of the period, with a value of 1.13 Jy 
on 2011 Feb 13 (day 926). During this period there was a fur- 
ther increase in the 15 GHz flux density at level of ~ 0.6 Jy. 
Instead, the activity at y-rays was sparse and of the order of 
a few times 10~ 7 ph cirT 2 s . NIR observations on April 26 
(day 998) indicated a high level (~ 1.2 mJy) that dropped by 
a factor 3 within a couple of weeks. About ten days earlier 
(day 989) there was some high-level activity at y-rays with 
fo.i-ioocev = (1.3 ± 0.7) x 10- 6 ph cm" 2 s" 1 . 

In the optical, the U flux density increased up to 0.3 mJy (day 
1037, Jun 4) and then decreased again. On Jun 13 (day 1046), 
the 15 GHz flux density reached the local maximum of ~ 0.7 Jy 
and a few days later, on Jun 20 (day 1053), the source erupted to 
another strong y-ray outburst (L y ~ 10 48 erg s _1 ). This event was 
reported by the AGILE satellite (AGILE Coll. 201 1). Another y- 
ray outbursts at a level of ~ 10~ 6 ph cirT 2 s _I occurred on Jul 28 
(day 1091). 

An interesting event occurred on October 9-12 (days 1164- 
1167). Swift observed the highest optical/X-ray fluxes ever 
recorded (Fo.3_ io keV ~ l.OxlO -11 ergcirT 2 s -1 ; F\j ~ 0.43 mJy), 
but with no corresponding high activity at y-rays ("orphan 
flare"). High flux density was recorded by IRAM at 86 GHz as 
well, while at lower frequencies (37 and 15 GHz) there was a 
decreasing trend. VLBA 2-cm observations at the end of 2011 
indicated a decrease of the polarisation and some hint for an- 
other swing in the EVPA. 

3.5. Time variability 

To obtain an estimate of the shortest observed variability in 
each band, we follow the method applied in, e.g., Foschini et 
al. (2011c). We calculated the observed time scale r for dou- 
bling/halving flux according to the formula: 

F(t) = F(t ) ■ 2- { '-' oVt (2) 

where F(t) and F(to) are the fluxes at the time t and fo, respec- 
tively, of two adjacent points. Although, it is more common to 
find in the literature the time scale for an exponential flux change 
(e.g. Scargle 198 1, LAT Coll. 2010e), we prefer to adopt the dou- 
bling/halving t because of the large error bars in the y-ray curve. 
The error in the evaluation of t has been calculated according to 
the bootstrap method, by assuming the worst case (T max ) calcu- 
lated by using the smallest flux difference and the longest time 
interval. The results are displayed in Table [2] and includes the 
shortest r and also the flux ratio measured between the maxi- 
mum and minimum of the selected curve. It is worth noting that 
the displayed r is the observed value. The intrinsic characteristic 
time scale can be obtained by dividing r by (1 + z). 

We note some peculiar variations in our estimates in Table|2] 
for example the wide variation between neighboring frequencies 
at 37 GHz and 32 GHz. While such a sensitivity of variability 
rates with frequency would be quite interesting, it is likely that 
this is an artifact of the unequal sampling. Effelsberg almost reg- 
ularly sampled the source with one-month cadence, while other 
curves, with more irregular sampling (e.g. 37 GHz), made it pos- 
sible to find also shorter variability (~ 2 days), depending on the 
clustering of observations around specific events. The most re- 
liable results are those at y-rays and at 15 GHz, which are also 
the curves most densely sampled. In these cases, it is possible to 
find variability on a few days time scale. The detection of such 
a day-scale variability at radio frequencies is intriguing, but not 
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Table 2. Variability at different wavelengths. Col. [1]: 
Band/Filter/Frequency of the data; Col. [2]: Observed charac- 
teristic time-scale for doubling/halving the flux [day] calculated 
over two consecutive points; Col. [3]: Significance of the flux 
variation [<x]; Col. [4]: Maximum flux ratio measured over the 
whole light curve. 
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59 ± 11 


8.8 


5.9 


86 GHz 


55 ± 14 


10.0 


4.0 


43 GHz 


236 ± 46 


13.1 


2.2 


37 GHz 


2.2 ± 1.8 


3.4 


4.3 


32 GHz 


83 ±23 


3.7 


4.4 


23 GHz 


< 244 


3.3 


2.5 


15 GHz 


< 6 


3.6 


5.8 


10.45 GHz 


83 ±6 


24.8 


3.5 


8.4 GHz 


40 ±22 


3.8 


3.1 


5 GHz 


117 ±42 


7.5 


3.6 


2.64 GHz 


179 ± 20 


20.0 


1.9 



extraordinary, since it was already detected in other blazars (cf. 
Wagner & Witzel 1995). 

The y-rays show the most extreme flux variations, with a 
time scale of a few days, although during outbursts - when the 
statistics are sufficient to perform an analysis with shorter time 
bins - it is possible to find intraday variability (Foschini et al. 
201 la). The relativistic jets of blazars are known to display sim- 
ilar time scales (cf. LAT Coll. 2010e), and during outbursts it is 
possible to measure variability on time scales of hours (Foschini 
et al. 2011c). The lowest flux changes are in the UV, which is 
expected to be dominated by radiation from the accretion disk 
as in the quasar subclass of blazars (cf Bonning et al. 2012), 
and at 2.64 GHz, which should refer to the largest spatial scales. 
Flux changes by factors 5-6 occur at X-rays, infrared H filter, 
142 GHz, and 15 GHz. To summarize, the observed variabil- 
ity (amplitude and time scale), by taking into account the bias 
of the sparse sampling, is consistent with that displayed by fiat- 
spectrum radio quasars. 

4. Spectral Energy Distribution (SED) Modelling 

We selected some specific cases to represent the behaviour of 
PMN J0948+0022, and for which we had sufficient data. It is 
worth noting that the availability of X-ray observations played an 
important selection effect: indeed, as shown in Fig. [4] the peak of 
the synchrotron emission cannot be determined accurately, be- 
cause of the synchrotron self-absorption. The X-ray emission is 
assumed to be due to the synchrotron self-Compton process; in 
this case, X-ray observations are required to estimate the mag- 
netic field strength for modelling purposes. Therefore, even if 
the light curves (Fig. [3]) indicate the presence of several interest- 
ing features, it was possible to properly model the corresponding 
SEDs only if an X-ray observation had been done within a few 
days. 



The selected time periods are (see Fig. [3]): 

- 2009 May 5: early decrease of the optical-to-X-ray flux dur- 
ing the 2009 MW Campaign (ferm//LAT data integrated 
over a time period of 5 days centred on the Swift observa- 
tion; data at other wavelengths at the closest time); 

- 2009 May 15: lowest X-ray flux recorded during the 2009 
MW Campaign and of the whole ~ 3.5 years (ferm//LAT 
data integrated over a time period of 5 days centred on the 
Swift observation; data at other wavelengths at the closest 
time); 

- 2010 July 8: first y-ray outburst (ferm//LAT data integrated 
over one day; Swift observation on July 3; data at other wave- 
lengths at the closest time); 

- 2011 June 21: second y-ray outburst (ferm//LAT data inte- 
grated over one day; Swift observation on June 14; data at 
other wavelengths at the closest time); 

- 2011 October 9-12: maximum optical-to-X-ray flux 
(ferm//LAT data integrated over a time period of 5 days 
centred on October 1 1 ; two Swift observations on October 9 
and 12 integrated; data at other wavelengths at the closest 
time). 

The corresponding SEDs and models are displayed in Fig. [4] 
It is worth noting that the data are not strictly simultaneous. As 
written above, the zero-order rule was to use the Swift observa- 
tion as the main driver, but obviously there were major events 
(e.g. the y-ray outbursts) that cannot be ignored. When possi- 
ble, we used the closest (within one week) Swift pointing, which 
were on 2010 July 3 for the 2010 outburst (i.e. 5 days before) 
and on 2011 June 14 for the 2011 outburst (i.e. 7 days before). 
In 201 1 there were other outbursts at y-rays (201 1 July 28, day 
1091; 201 1 September 25, day 1 150), but no timely Swift obser- 
vations were available. 

We adopted the same model used in past MW campaigns 
(LAT Coll. 2009c, Foschini et al. 2011), which is that devel- 
oped by Ghisellini & Tavecchio (2009). It is a one-zone model, 
which calculates the emitted radiation from a population of rel- 
ativistic electrons through synchrotron self-Compton and exter- 
nal Compton (EC) processes. In the latter, the seed photons are 
considered from different sources (accretion disk, broad-line re- 
gion, infrared torus, etc.). We refer to Ghisellini & Tavecchio 
(2009) for more details on the model and to LAT Coll. (2009c) 
and Foschini et al. (201 1) for the specific application to the 2009 
and 2010 MW campaigns on PMN J0948+0022, respectively, 
with specific reference to the possible changes in the model pa- 
rameters as a function of the data. 

The input parameters (see Table [3]) were adjusted to fit the 
data and the model output consists of the calculated jet power 
divided into its basic components (radiative and kinetic, see 
Table 0J. Given the large number of parameters, we fixed some 
of them by taking into account the available measured quanti- 
ties. The mass of the central black hole is M ~ 1.5 x 1O 8 M , 
as derived from the fit to the accretion disk emission (LAT Coll. 
2009b), which is consistent with the black hole mass estimate 
by Zhou et al. (2003) by using the classical virial method. The 
accretion disk luminosity is measured as Ld ~ 9 x 10 45 erg s , 
corresponding to Ld/^Edd ~ 0.4. The value has been measured 
during the period of minimum jet flux and fixed on short time 
scales, although it can change over the years. 

The jet viewing angle is fixed to ~ 3° for all SEDs. In earlier 
work, we adopted a value of 6° (LAT Coll. 2009b, c), but mod- 
elling of the 2010 July outburst required a smaller viewing an- 
gle to explain the observed y-ray isotropic luminosity (Foschini 
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Fig. 4. {left panel) Spectral Energy Distributions (SEDs) of PMN J0948+0022 in some specific time periods, (right panel) Zoom of 
the infrared-to-y-ray energy range. Each period is indicated with different symbols and colours, as displayed in the figures. Since 
we adopted an one-zone model, radio data were not fitted. The bow-ties at y-rays indicate only the changes in the photon index, 
which are most relevant for the modelling. It is necessary to take into account that the flux normalisation has an error too, although 
not displayed. See the text for details. 



Table 3. Input parameters used to model the SED. Col. [1]: Date. Col. [2]: distance of the blob from the black hole in units of 10 
cm; Col. [3]: /?diss in units of the Schwarzschild radius Rs ; Col. [4]: power injected in the blob calculated in the comoving frame, in 
units of 10 45 erg s _1 ; Col. [5]: magnetic field in Gauss; Col. [6]: bulk Lorentz factor; Col. [7]: random Lorentz factor of the injected 
electrons at the break of the distribution; Col. [8]: maximum random Lorentz factor of the injected electrons; Col. [9] and [10]: 
slopes of the injected electron distribution below and above yt,. All models assume a black hole mass M = 1.5 X 1O 8 M and 

an accretion disk luminosity Ld = 9 x 10 45 erg s _1 , corresponding to La/Lndd = 0.4. The X-ray corona has a luminosity Lx = 0.3Ld- 
The radius of the BLR is fixed to Rblr — 3 x 10 17 cm. The spectral shape of the corona is assumed to be oc v exp(-/iv/150 keV). 
The viewing angle is 3° for all models. 
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13 


10 
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2.65 


2009 May 15 


49.5 


1100 


0.0055 


6.0 


11 


50 


1500 





2.7 


2010 Jul 08 


103.5 


2300 


0.044 


2.1 


16 


50 


3000 





2.8 


2011 Jun21 


108.0 


2400 


0.025 


2.2 


16 


240 


4000 


0.5 


3.9 


2011 Oct 09 


90.0 


2000 


0.032 


3.3 


12 


80 


5000 


1 


3.05 



et al. 2011). The smaller jet angle generates slightly different 
parameter values compared to LAT Coll. (2009c), but does not 
significantly change the interpretation. 

In addition, fixing the above parameters has some impli- 
cations on the other parameters. For example, as known from 
the reverberation mapping, the radius of the broad-line region 
(BLR) is a function of the accretion disk luminosity (see Eq. 2 
in Ghisellini & Tavecchio 2009). The values of the fitted parame- 
ters are presented in Tableland the calculated output jet powers 
are shown in Table [4] 

The modelled SEDs can be interpreted by looking at the most 
relevant changes in the parameters. The drop in y-ray emission 
in early 2009 May can be basically explained with a decrease in 
the injected power (from 3 x 10 43 erg s _I to 5.5 x 10 42 erg s~'), 
an increase of the magnetic field (from 4.4 to 6.0 G), and a de- 
crease of the bulk Lorentz factor (from 13 to 11). The slope 
of the electron distribution remained almost constant. On the 
other hand, the 2010 July outburst requires a high injected power 
(4.4 x 10 43 erg s~') and also a higher bulk Lorentz factor (F = 
16), together with a smaller magnetic field (2.1 G). The 2011 



Table 4. Calculated jet power in the form of radiation, Poynting 
flux, bulk motion of electrons and protons (assuming one cold 
proton per emitting electron). The logarithms are calculated on 
the powers expressed in erg s _1 . 



Date 


log/> r 




logPe 


log^p 


2009 May 05 


45.31 


44.75 


45.23 


47.49 


2009 May 15 


44.55 


44.59 


44.33 


46.17 


2010 Jul 08 


45.97 


44.67 


45.43 


47.61 


2011 Jun21 


45.76 


44.47 


44.98 


47.10 


2011 Oct 09 


45.47 


44.68 


45.20 


47.46 



June outburst, although with similar magnetic field and Lorentz 
factor (2.2 G and 16, respectively), requires less injected power 
(2.5 x 10 43 erg s~'), but a softer electron distribution. Indeed, as 
noted in Sect. 2.1, the LAT spectrum showed in 2011 softening 
in the slope. 

A comparison of the values obtained by the SEDs model- 
ing of PMN J0948+0022 (Table [3] and © with those of a large 
sample of bright blazars (Ghisellini et al. 2010) shows that the 
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jet powers are in the range of flat-spectrum radio quasars, al- 
though there are some differences in the input parameters. The 
average mass and accretion luminosity of quasars in the sam- 
ple of Ghisellini et al. (2010) are ~ 1O 9 M and ~ 0.1L Edd , re- 
spectively, to be compared with the values of PMN J0948+0022 
(M ~ 1O 8 M , L d ~ 0.4L Edd ). However, it is possible to find 
quasars with similar characteristics, like PKS 0426-380 (z = 
1.112), which has M ~ 4 x 10 8 M o and L d ~ 0.6L Edd . The in- 
put and output parameters of the model (see Table 4 and 5 in 
Ghisellini et al. 2010) are very similar to those presented here. 

5. Conclusions 

The most likely physical interpretation for the observed data is 
that the NLS1 galaxy PMN J0948+0022 hosts a powerful rela- 
tivistic jet directed, within a few degrees, toward the Earth, con- 
firming previous works. Zhou et al. (2003) and Doi et al. (2006), 
having found an inverted radio spectrum with high brightness 
temperature in excess of 10 13 K, suggested that the NLS1 galaxy 
should host a relativistic jet. However, there are some caveats 
in inferring the presence of a relativistic jet from the high- 
brightness temperature, as outlined for example by Tsang & Kirk 
(2007) and Singal (2009). A halo of doubt remained, particu- 
larly given the fact that it was quite anomalous and strange that 
a NLS 1 type AGN, which is generally radio-quiet, hosted a rel- 
ativistic jet. 

The detection of variable high-energy y-ray emission from 
PMN J0948+0022 with Fermi/LAT (LAT Coll. 2009b, 2010a) 
provided the needed breakthrough. Specifically, the 2009 MW 
campaign, which found coordinated MW variability, confirmed 
beyond any doubt the association of the high-energy y-ray 
source with the NLS1 (LAT Coll. 2009c). In addition, the ob- 
servation of optical (V filter) polarisation at 19% level (Ikejiri et 
al. 201 1) and optical (B and R filters) intraday variability (Liu et 
al. 2010), provided useful and important complementary infor- 
mation to strengthen the relativistic jet scenario. 

The reanalysis of the data presented in this work confirm 
and extend the early findings reported by LAT Coll. (2009b, c, 
2010a). We have selected five epochs to represent the source ac- 
tivity during these three and a half years. The most relevant pa- 
rameters driving the modelling of the SEDs in 2009 and 2010 
are the injected power, the bulk Lorentz factor, and the mag- 
netic field. The 201 1 SEDs are instead characterised by a change 
(softening) in the electron distribution, with a decreasing contri- 
bution of high-energy leptons. The analysis of the SEDs mod- 
elling, together with the study of the variability, indicates that the 
jet hosted by the NLS 1 PMN J0948+0022 is like those hosted by 
the flat-spectrum radio quasar subclass of blazars. This is another 
point favouring the universality of the jet phenomenon. 

In the years 2009 and 2010, the jet emission followed what is 
expected from the canonical model of relativistic jets (Blandford 
& Konigl 1979, Heinz & Sunyaev 2003), where the peak of ra- 
dio emission follows the y-rays after a few months. The year 
2011 was instead characterised by a prolonged activity, with 
some peculiar events. Specifically, there was an optical/X-ray 
flare, which has no corresponding activity at other wavelengths. 
Such "orphan flares'0 have been observed also in other blazars 
(e.g. Krawczynski et al. 2004, LAT Coll. 2010d, VERITAS 

8 Although this term has been created by H. Krawczynski when 
speaking about a VHE flare of 1ES 1959+650 with no counterpart at 
other wavelengths (Krawczynski et al. 2004), its meaning can be ex- 
tended to any flare at any frequency with no counterparts at other fre- 
quencies. 



Coll. 2011, Marscher 2012, Marscher et al. 2012). Some at- 
tempts to explain this behaviour have been suggested (e.g. 
Marscher 2012), but none seems to meet with the general con- 
sensus. In the present case, the sparse and irregular sampling at 
NIR/optical/UV/X-rays prevents the testing of specific hypothe- 
ses. The observed variability of a few days suggests the need of a 
denser sampling rate, perhaps a MW campaign limited to a short 
period, but with day-scale sampling. 
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Appendix A: Light curves 

Light curves of all the data available and studied in the present 
work are displayed in the following. 
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Fig. A.l. Light curves at various frequencies. From top left to bottom right panel: y-ray 0.1 - 300 Ge V from Fermi/LA^T, with 1 day 
time bin [10~ 7 ph cirT 2 s -1 ] (the horizontal dotted line correspond to a TS =4 average upper limit of 2.3 x 10~ 7 ph crrT 2 s _1 for a 
one-day exposure); X-ray 0.3 - 10 keV from Swift/KKT [c s 1 ]; UVW2 from Swift/UVOT [mJy]; UVM2 from Swift/UVOT [mJy]; 
UVW1 from Swift/UVOT [mJy]; U from Swift/UVOT [mJy]. Time starts on 2008 August 1 00:00 UTC (MJD 54679). 
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Fig. A.2. Light curves at various frequencies. From top left to bottom right panel: B from Swift/UVCfT [mJy]; V from Swift/UWOT 
[mJy]; J from INAOE [mJy]; H from INAOE [mJy]; K from INAOE [mJy]; 142 GHz from IRAM [Jy]; 86 GHz from IRAM [Jy]; 
43 GHz from Effelsberg [Jy]. Time starts on 2008 August 1 00:00 UTC (MJD 54679). 
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Fig.A.3. Light curves at various frequencies. From top left to bottom right panel: 37 GHz from Metsahovi [Jy]; 32 GHz from 
Effelsberg [Jy]; 23 GHz from Effelsberg [Jy]; 15 GHz from OVRO (open squares), Effelsberg (filled circles), and MOJAVE (open 
triangles) [Jy]; 10.45 GHz from Effelsberg [Jy]; 8.4 GHz from Effelsberg (filled circles) and Medicina (filled stars) [Jy]; 5.0 GHz 
from Effelsberg (filled circles) and Medicina (filled stars) [Jy]; 2.64 GHz from Effelsberg [Jy]. Time starts on 2008 August 1 00:00 
UTC (MJD 54679). 



